Breaking and restoration of rotational symmetry for irreducible tensor operators on the lattice 
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We study the breaking of rotational symmetry on the lattice for irreducible tensor operators and practical 
methods for suppressing this breaking. We illustrate the features of the general problem using an a cluster model 
for ^Be. We focus on the lowest states with non-zero angular momentum and examine the matrix elements 
of multipole moment operators. We show that the physical reduced matrix element is well reproduced by 
averaging over all possible orientations of the quantum state, and this is expressed as a sum of matrix elements 
weighted by the corresponding Clebsch-Gordan coefficients. For our a cluster model we find that the effects of 
rotational symmetry breaking can be largely eliminated for lattice spacings of n < 1.7 fm, and we expect similar 
improvement for actual lattice Monte Carlo calculations. 

PACS numbers: 12.38.Gc, 03.65.Ge, 21.10.Dr 




I. INTRODUCTION 


In recent years, lattice Monte Carlo calculations have been widely 
■^pplied to the study of nuclear structure In particular, chiral 

^ffective field theory (chiral EFT) combined with lattice methods 
^as been employed to study the spectrum and structure of light and 
'—irledium-mass nuclei [Q-^. In such calculations, continuous space- 
t^me is discretized and compactified so that path integrals can be 
^omputed numerically. The mesh points uniformly span a cubic 
ly-^ox, and some boundary conditions such as periodic boundaries are 
O^nposed in each dimension. However, the calculated bound state 
kOr.ergies and wave functions will, in general, deviate from their con- 
^TSnuum infinite-volume counterparts due to the discretization and 
"^^nite-volume artifacts. 

Over the years, much effort has been devoted to removing nu- 
•^^erical artifacts in lattice field theory calculations. The finite vol- 
lOrne energy shifts for two-body bound states [^ -ITsl ] as well as two- 
Ttody resonances and scattering problems [[l9|^3[] have been stud- 
^ied in detail. There is also on-going research to extend these results 
kT( 3 bound states with more than two constituents [ p^ [T3| , |24|-|^. 
the other hand, removing the artifacts from finite lattice spacing 
a more complicated issue. For chiral EFT the lattice improve¬ 
ment program proposed by Symanzik et al. [|2^-|^ provides a use¬ 
ful approach for systematically reducing discretization errors. This 
method was also applied to Yang-Mills theories gauge field 

theories [^l|-p^ and QCD [^Sj]. Dudek et a l. have proposed a 
method where the continuum spin of meson []37[] and baryon [ |^ , 
excited states in lattice QCD can be reliably identified. Meanwhile, 
Davoudi et al. have quantified the breaking and restoration of 
rotational invariance at both tree level and one-loop level by means 
of lattice operators smeared over a finite spatial region. 

On the lattice, the full rotational symmetry group is reduced to 
the finite group of cubic rotations, SO{3) -A 5(9(3,Z). Several ba¬ 
sic rules based on the argument of rotational invariance are broken. 
For example, in the continuum and infinite-volume limits, quantum 
bound states with angular momentum J form a degenerate multiplet 
consisting of 27 -F 1 components, while on the lattice the energy lev¬ 
els split into subgroups corresponding to different irreducible repre¬ 
sentations {irreps) of the cubic group iMl- The size of the en¬ 


ergy splittings are dictated by the lattice spacing and by the volume 
and boundary conditions. 

In Ref. iQ], we explored the breaking of rotational symmetry 
on the lattice for bound state energies with an a cluster model. It 
was shown that the calculated energy is minimized when the natural 
separation between particles is commensurate with the separation 
between lattice points along the preferred lattice directions associ¬ 
ated with the given angular momentum state. It was also shown that 
the multiplet-averaged energy is closer to the continuum limit than 
any single energy level. One can apply these results to future ab ini¬ 
tio lattice simulations for nuclear systems where a cluster structures 
are important. 


In this paper, we extend the analysis of Ref. |44| to other ob¬ 
servables besides the energy. One is often interested in the nuclear 
radii, quadrupole moments as well as transition probabilities. For 
example, an anomalously large radius compared with the usual A* 
scaling law is evidence for a halo nucleus while the in¬ 

trinsic quadrupole moment is often related to the rotational bands 
observed in deformed nuclei [BS]. 


We consider irreducible tensor operators sandwiched by a pair of 
bound state wave functions. In the continuum limit such an expres¬ 
sion can be factorized and simplified according to the Wigner-Eckart 
theorem because of the full rotational symmetry, but on the lattice 
such factorization is no longer possible and the situation becomes 
more complicated. Furthermore, continuum selection rules for elec¬ 
tromagnetic transitions are not exactly satisfied on the lattice. Some 
transitions that are absolutely forbidden by rotational symmetry may 
assume non-zero amplitudes on the lattice. We would like to con¬ 
struct proper corrections for these matrix elements in order to mini¬ 
mize symmetry breaking effects. 


Our objective in this paper is to investigate anisotropic lattice ar¬ 
tifacts in the matrix elements of irreducible tensor operators and 
search for a practical method to restore rotational symmetry. The 
details of the particular interaction are not essential to our general 
analysis. Here, we use the same a cluster model as in Ref. [Q] 
where the a-a interaction is approximated by an Ali-Bodmer type 
potential adjusted to produce a bound ^Be nucleus. 






2 


II. THEORETICAL ERAMEWORK 


B. Lattice wave functions 


A. Hamiltonian 

Let ma = 3727.0 MeV denote the mass of the a particle and 
m = ma/2 the reduced mass. Our starting point is the one-body 
Hamiltonian 

H=-|^ + VM, (1) 

where r = is the distance between the two a particles, 

V = VV -I- Vc is the a-a potential, including nuclear and Coulomb 
potentials. 

We use the same two-body potential as that used in Ref. 0. 
For completeness we briefly introduce its functional form and the 
parameters. For the nuclear part of the a-a interaction we use an 
isotropic Ali-Bodmer-type potential. 


VN{r) =yoexp(-77or2)-)-yiexp(-77fr2), (2) 

where Vq = —216.0 MeV, Vi = 354.0 MeV, rjo = 0.436 fm^' and 
rji = 0.529 fm^*. These parameters are determined by fixing the 
S- and D-wave a-a scattering lengths to their experimental values. 
The Coulomb potential is given by the error function 


4£?^ 

Vc(r) = —erf 



( 3 ) 


where Ra = 1.44 fm is the radius of the a particle, e is the funda¬ 
mental unit of charge and erf denotes the error function. 

Since our focus is on measuring bound state properties and the 
physical ^Be nucleus is unbound, we increase Vq by an amount 
of 30%. With this strengthened potential, the ^Be nucleus has 
a ground state at £(0+) = —10.8 MeV and one excited state at 
£(2+) = —3.3 MeV. These energies are measured relative to the 
a-a threshold. 

In the lattice calculations the spatial vector r assumes discrete 
values, and so the Hamiltonian in Eq. (1) becomes a matrix. Pe¬ 
riodic boundary conditions for a box of size £ are imposed on the 
wave functions. 


V/(/- + n,£) = v/(r), (4) 

where n,- with i = x,y,z are the unit vectors along the three coordi¬ 
nate axes. The energy eigenvalues and wave functions can be ob¬ 
tained by the diagonalization of the Hamiltonian matrix of dimen¬ 
sion 

The kinetic energy term in Eq. (1) can be expressed on the lattice 
by finite differences. In one dimension we have 

/"W [/(x-fka)-f/(x-ka)], (5) 

k=l 

where a is the lattice spacing and is a set of fixed coefficients. 
The order-V formula involves 2N 1 points and the corresponding 
truncation error is In Ref. 0], we gave the coefficients 

4^' forV < 4. In this paper, we use the V = 4 formula to calculate 
the second-order derivatives appearing in the Laplace operator. This 
choice removes most of the rotational symmetry breaking effects 
due to the difference formula for the kinetic energy. 


The continuum Hamiltonian Eq. (|I|) is invariant under any spatial 
rotation. As a result, the bound states of H form angular momen¬ 
tum multiplets. Let us denote the bound state wave functions as 0/„,, 
where the integer I is the angular momentum and the integer m is 
the z-component of angular momentum with —I < m < 1. Eor sys¬ 
tems with more than one bound state with the same value of I, we 
will need additional radial quantum numbers. However, we do not 
consider such cases here. The angular dependence of these wave 
functions are given by the corresponding spherical harmonics Yi„. 

On the lattice, these angular momentum I multiplets are split into 
irreps of the cubic rotational group, 50(3,Z). In Ref. [0, we gave 
the splitting patterns of the angular momentum multiplets with / < 8. 
In order to specify the wave functions belonging to the same irrep, 
we define a quantum number k that is valid on the lattice through 
the relation 


K.(f)=exp(-i|t), (6) 

where R^{n/2) is a rotation around the z-axis by ;r/2. The integers 
k are equal to m modulo 4 and are non-degenerate for each irrep of 
O. We label the wave function for any eigenstate by /, k and 
the irrep T it belongs to. If the angular momentum J contains more 
than one branch belonging to the same irrep, we distinguish them by 
adding primes to the name of the irreps. Eor instance, the symbol 
denotes the wave function with I = 6, k = 1 belonging to 
the second T 2 irrep. Consequently, in the continuum limit the wave 
functions form a complete basis for the bound state subspace, 
and the corresponding energies are degenerate for the same total 
angular momentum J, while on the lattice the energy eigenvalue 
depends on both J and the irrep z. In the continuum limit we can 
write down unitary transformations from the basis to the vz/tA: 
basis and vice versa, 

0/m — ^^^ImxkWlTk: (7^ 

tk 

WHk = Y.^HknAn- (8) 

m 

See Ref. [ 0 ] for details of these transformations. As an example, 
we show the case of / = 2. The wave functions V^ 2 £ 0 j W 2 E 2 belong 
to irrep £, and 1 / 2721 ’ ¥ 2 t, 2 , W 2 T )3 belong to irrep T 2 . Eollowing 
Ref. [0 we find 

¥2E0 = 020, W 2 T 2 I = 021, W 2 T 23 = 02l, 

¥2E2 = Y^(022 + 022) V 2 T 22 = -*^^(022 “ 022)' (9) 

We write m for notational convenience to denote —m. 

On the lattice we can obtain the bound state wave functions 
by simultaneously diagonalizing the lattice Hamiltonian H (or trans¬ 
fer matrix [|I]]) and the Rt{7z/2) operator. Since the full rotational 
symmetry is broken, the angular momentum / should be viewed as 
a label that describes the angular momentum multiplet obtained by 
dialing the lattice spacing continuously to zero. But we can use the 
unitary transformation in Eq. to define the wavefunctions 0;^ at 
non-zero lattice spacing. We do this even though the wavefunctions 
0/,„ are generally not exact eigenstates of H when the lattice spacing 
is non-zero. 
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Consider now the bound state wavefunctions for a zero angular 
momentum state 0oo and a general state 0/,„. In the continuum limit 
the matrix elements of a r^Yi,„ multipole operator inserted between 
000 and must be independent of m, 

{0lm 1^ Yjffj 10oo) 

= J c/V(0;,„|p(r)/y,,„(n)|0oo) =c„, -;>C. (10) 

When using Eq. (^, we find the condition that C,„ is independent 
of m a very convenient check that the phases for 0/„, and are 
consistent with standard conventions as defined in Ref. [p3|]. 


C. Factorization of the matrix elements 


Given a pair of bound state wave functions (r) and (r) 
with angular momenta li, nii, I 2 and m 2 on the lattice, we can eval¬ 
uate the matrix element of the multipole moment operator 

{l\m\\r 0/jm| ^/m(tt) 0 / 2 m 2 (11) 

n 

where n runs over all lattice sites. Here, we choose independent 
integers / and /' in the multipole moment operator, in order to keep 
the radial and angular degrees of freedom independent. This makes 
our conclusions sufficiently general and applicable to all irreducible 
tensor operators. We use parentheses to denote matrix elements on 
the lattice, {f\0\i), and Dirac brackets to denote matrix elements in 
the continuum limit, {f\0\i). The continuum limit of Eq. dill ) is 

{hmi\/Yi„\l2m2) = J iiV0,*,„j(/-)/T/,„(n)0,2m2('') (12) 

where the integration is performed over all space. Matrix elements 
in the form of Eq. ( [I^ occur frequently in the calculation of vari¬ 
ous nuclear observables such as the mean square radii, quadrupole 
moments, transition probabilities, etc. Here we focus on lattice ar¬ 
tifacts that produce some residual difference between the values in 
Eq. (|i^ and Eq. (|^ and the methods for removing them. 

Let us start with the continuum limit. According to the Wigner- 
Eckart theorem, the matrix element ( 0 ) can be transformed into 
a product of two factors: Clebsch-Gordan (C-G) coefficients and 
the reduced matrix element that contains the essential non-trivial 
physics. Let Ri{r) and R 2 ir) be the radial parts of the wave func¬ 
tions (r) and 0 /,m 2 (^)’ respectively, so that 

(Iimi|/y,„,|l2m2) = {h\/\l2)Q',lZli,n, (13) 

where 

{iy\l2)= J dr/+^RUr)R2{r) (14) 

QiZlln. = I (15) 

The radial integral in Eq. ( p^ is the matrix element of the /'-order 
moment operator and is independent of the quantum numbers nii, 
m and m 2 . Meanwhile can be written as a product of C-G 

coefficients. 


Qhm 

^/2m2,//n 


(2/ -I- 1)(2/2-|- 1) ;,o /|,„j 

47 r(2 /i-fl) bO,/o hnnM- 


( 16 ) 


All of the dependence on the quantum numbers mi, m, m 2 and / is 
absorbed into Instead of simply the factoring out the C- 

G coefficient cj";> we find this factorization useful because the 


radial integrals of r^' are physically meaningful. In Table |[ we list 


some factors Q 


/i»ii 
l2tn2,lm 


with 1 1 =12=2 and / from 0 to 4. The other 


factors can be obtained from standard tables of C-G coefficients. 

Since full rotational symmetry is broken on the lattice, it is not 
possible to write the wave functions as products of radial and angu¬ 
lar parts. However, one can still apply the Wigner-Eckart theorem to 
each irrep of the cubic group. As a result, the matrix elements calcu¬ 
lated on the lattice belonging to the same irreps are related by C-G 
coefficients of the cubic group, which can be computed easily using 
decompositions into spherical harmonics ^|. Eor example, we 
have 


(272l|r2y2£o|272l) = -y|(27’2l |r2y2£2|27’2l), (17) 

where y 2 £o Y 2 E 2 are defined as 


Y 2 EQ = Y 20 , Y 2 E 2 = \j^iY22+Y22), (18) 

in analogy to the relations in Eq. (^). The ratio of — a/ 1/3 is irre¬ 
spective of the lattice spacing, box size or the strength of the inter¬ 
action. 

We will divide the lattice matrix elements in Eq. ( Ill] ) by the 
Qk^im as defined i^q. ( p^ whenever non-zero, even though the 
factorization in Eq. ( jT^ is not exact on the lattice. We represent the 
resulting quantity by double vertical lines, 

(/imi||/y;,„||/2m2) = ihmi\/Yi„,\l2m2)/Q‘,'ZZZim^ (19) 

for 2 / 2 m 2 im 7^ clear that these reduced matrix elements all 

converge to the radial matrix element (/i \r''\l 2 ) as a —>■ 0. At non¬ 
zero lattice spacing, however, the ratio will depend on the quantum 
numbers mi, m 2 and m. Therefore, the splittings between differ¬ 
ent components of Eq. (|^ are sensitive indicators of the rotational 
symmetry breaking effects. 


D. Isotropic average 

We now focus on spatial anisotropies which are associated with 
the orientation of our lattice wave functions 0 /,„ relative to the lat¬ 
tice axes. Let us illustrate this point by recalculating the matrix el¬ 
ement (/imi|r^ T/„, |/ 2 m 2 ) with a tilted lattice which differs from the 
original one by a rigid body rotation. Let us assume that the wave 
functions and 0 / 2 m 2 can be smoothly interpolated in between 
lattice points as a function in continuous space. Then we can define 
a tilted matrix element as 

(/imi|r' y/,„|/2m2)A 

n 

where A = (a, j3, 7 ) is a set of Euler angles and R{A) is an element 
of the SO(3) rotation group. In the continuum limit, rotational in¬ 
variance guarantees that (/imi |/ y/m|^ 2 '« 2 ) and (/imi \r^ y)m|^ 2 '« 2 )A 
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h I h 
2 0 2 


2 2 2 


2 4 2 


m\ m m2 
0 0 0 
2 0 2 

0 0 0 
2 0 2 
2 2 0 

1 2 T 

0 0 0 

2 0 2 
2 1 1 
1 1 0 
2 4 2 


Wl'Wl 

1 

IsfH 



mi m m2 
1 0 1 

2 1 1 
1 1 0 
1 0 1 

1 3 2 

1 2 T 
1 0 1 

2 2 0 


Qhm 

1 

Isfii 


/ 15 

Y 14;t 
1 /J 
T^Y ^ 

X /i 
14 Y 


i /x 

i\j 111 
_\ pM 

7\ JT 

_2_ 

iVn 

1 /TJ 
"H Y ^ 


Table I. The factor defined in Eq. (|i^) with li =12 = 2 and / from 

0 to 4. The other ratios not listed here can be obtained from standard tables 
of C-G coefficients. 


are equal. For non-zero lattice spacing and non-zero angular mo¬ 
menta, however, the tilted matrix element will depend on the quan¬ 
tity A. We would like to eliminate this unphysical orientation de¬ 
pendence from the final results. 

A natural choice is to average the results over all possible ori¬ 
entations, or equivalently, over the whole SO(3) group space. The 
isotropically-averaged matrix element is defined as 


{hmi\/Yi„\l2m2)o = j cP K{limi\/Yi,n\l2m2)K 




— c 


Jllfll 

l2m2,lin 


2/1 + 1 f , , ‘-"'2 

m' .m j ,m2 


. ( 21 ) 


Here, pA is the normalized invariant measure on the SO(3) group 
space and the subscript “o” denotes the isotropic average. In anal¬ 
ogy with Eq. ([T^, it is convenient to define a reduced isotropically- 
averaged matrix element, 

ihWr’YiWkP = ( 22 ) 


for im ^ continuum limit (/i ||/T/||/ 2 )o coincides with 

the radial matrix element (/i |l 2 )- 

We note that the radial matrix element (/i |r^ I/ 2 ) is not only inde¬ 
pendent of m, mi and m 2 , but also independent of 1. So a non-trivial 
test of rotational symmetry restoration is to check that {li\\P T;||l 2 )o 
as defined in Eq. ( p^ is independent of 1. If {h\\r' T/||/ 2 )o is to a 
good approximation independent of I, then we have succeeded in 
approximately factorizing the radial and angular parts of the lattice 
wave function by means of isotropic averaging. We will test this 
numerically with the a cluster model in Sec. |i|. 


m. RESULTS 


We start with the mean square radius operator P. This corre¬ 
sponds to setting I = 0 and /' = 2 in Eq. (|rT]) and Eq. ([1^. In the 
upper panel of Eig. we show the expectation values of P for the 
lowest 2+ states as functions of the lattice spacing a. The eigenstate 
wave functions 1 / 2 ^ are obtained by simultaneously diagonalizing 
the Hamiltonian H in Eq. (|^ and the R^k/I) operator. We then 
construct the linear combinations 02m according to Eq. (|9|). We write 
(m||0||m) as an abbreviation for (2m||r^yool|2nj) and the solid curve 
denotes the isotropic average (2||r^yool|2)o defined inEq. (p^. Only 
three values with m > 0 are shown since time reversal symmetry en¬ 
sures equal results for m and —m. As discussed in Sec. y, all these 
reduced matrix elements converge to {P) = {h = '^\P\h = 'P) as 
a ^ Q. Note that in all the following calculations, we remove the 
finite volume effects by using a box size of L > 16 fm. 

The three branches are not linearly independent because of the 
cubic symmetries on the lattice. Eor scalar operators, the linear re¬ 
lations among them are manifest. According to Eq. (^, the wave 
functions 021 and 02o belong to the irreps E and T 2 , respectively. 
Meanwhile the wave function 022 is an mixture of the irreps E and 
T 2 with equal weights. Thus, (2||0||2) equals the arithmetic average 
of (0||0||0) from irrep E and (1||0||1) from irrep 72- So we find 
that the isotropically-averaged reduced matrix element (2||r^yo||2)o 
is given by 

(2||r2yo||2)o = ^(l||0||l) + ^(0||0||0), (23) 


where the numerators 3 and 2 are simply the dimensions of each 
cubic representation. It is easy to verify that the weighted average 
formula is applicable for any angular momentum with the factors 3 
and 2 in Eq. ( |2^ replaced by the corresponding irrep dimensions. 
In Ref. iQ] we introduced the multiplet-weighted average to elim¬ 
inate the anisotropic effects in the bound state energies. We have 
now proven that this procedure is equivalent to averaging over lat¬ 
tice orientations and applies to the expectation value of any scalar 
operator. 

Next, let us examine the dependence of these reduced matrix el¬ 
ements on the lattice spacing a. Eor a < 1.0 fm, the three branches 
merge, and for large a they split and show oscillations. Before dis¬ 
cussing the details and the physics behind them, it is interesting to 
compare Eig. |I] in this paper to Eig. 3 in Ref. where the cal¬ 
culated 2+ energies are shown as functions of a. We immediately 
find that these figures look similar if we map (11|0|| 1) to £(2^), and 
(0||0||0) to £( 22 ). More specifically, in Eig. [| the splitting between 
the two branches has two zeros near a = 1.4 fm and 1.9 fm. Eor 
a < 1.0 fm, the splitting is negligible. Eor 1.0 < a < 1.4 fm, (1110111) 
is higher than (0||0||0). However, in the region 1.4 < a < 1.9fm, the 
order is reversed. Eor a > 1.9 fm, (1110111) is once again higher and 
the splitting increases monotonically. This behavior also occurs for 
the energies in Ref. |44| with slightly different turning points. Ad¬ 
ditionally, the weighted averages (2||r^yo||2)o and£(2^) both show 
down bending in the transitional region 1.5 < a < 2.0 fm, result¬ 
ing in smaller expectation values for both energy and radii at large 
lattice spacings. 

The same pattern is also observed for other scalar operators. In 
the lower panel of Eig. [|, we show the results for the r^ operator 
which corresponds to setting I = 0 and /' = 4. Here, the symbols 
(ot||0||ot) are abbreviations of (2OT||r^yoo||2OT) and the solid curve 
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Figure 1. (color online). Upper panel: Mean square radii (r^) for the lowest 
2+ states of the *Be nucleus as a function of a. (a||0||o!) is an abbreviation 
of the reduced lattice matrix element {2a\\r^YQQ\\2a) defined in Eq. (j^, 
which converges to (r“) as a ^ 0. The box size L is kept larger than 16 
fm to remove finite volume effects. The solid line represents the isotropic 
average (2||r^To||2)o defined in Eq. Lower panel: Mean value (r^) for 
the lowest 2+ states of the ^Be nucleus as a function of a. (a||0||a) is an 
abbreviation of the reduced lattice matrix element (2a||r^Tool|2«), which 
converges to (r"*) as a —> 0. The solid line represents the isotropic average 

(2|k4yol|2)o. 


denotes the isotropic average. All curves converge to the expecta¬ 
tion value (r'*) = (/i = 2 |r^|l 2 = 2) as 0. Again, the isotropic av¬ 
erage is just the multiplet-weighted average over the five-fold mul- 
tiplet 02 h;- Apparently, the oscillations of these components as well 
as the down bending of the isotropic average are quite similar to the 
ones observed for the operator. 

Based on the results shown in Fig. we can compare the various 
components (ot|| 0 ||ot) to the continuum limit and select the ones 
with the least lattice spacing dependence. It seems that (1||0||1) 
and ( 0 || 0 || 0 ) are not optimal for both and r"* operators because of 
their large deviations from the continuum values for sparse lattices. 
The arithmetical average (2||0||2) and the multiplet-weighted aver¬ 
age are close to each other for all lattice spacings considered here. 
In principle, both of them can be used as a good approximation to 
the continuum limit. Nevertheless, the multiplet-weighted average 
is theoretically preferable because of its clear physical interpretation 
as isotropic averaging. 

Our conclusions for the operators r^, r'^ as well as the energy 
can be generalized straightforwardly to other scalar operators on 
the lattice. For example, the linear relation among the components 


(0||0||0), (11|0|| 1) and (2||0||2) is also satisfied there. In order to es¬ 
timate the continuum value, we can calculate the isotropic average 
according to Eq. which equals the weighted average over the 
angular momentum multiplet in the case of scalar operators. 

We now turn to the case of 1 = 2. In the upper panel of Fig. |[ the 
symbols (a|| j3 1 | 7 ) are abbreviations of the reduced matrix elements 
( 2 a||r^} 2 j 3 l| 27 ) defined in Eq. (p^, where the subscripts a, j5 and 
7 run from —2 to 2 and only the components with a = j3 -f 7 are 
shown. The solid curve represents the isotropic average defined in 
Eq. (p^. All these curves converge to the mean square radius (r^) 
as fl —0. Note that while the 27-1-1 wave functions in an angular 
momentum multiplet are mixed with each other on the lattice, some 
components with + y survive for large lattice spacings. How¬ 

ever, because the corresponding C-G coefficients vanish in this case, 
these components do not contribute to the isotropic average. 

Compared to the case of the scalar operator shown in Eig. [I|, 
the insertion of the spherical harmonic I 2/3 makes Eig. ^much more 
complicated. Nevertheless, we can still draw some general conclu¬ 
sions. Eirst, as for the scalar operators, we can show that (2||0||2) is 
just the arithmetical average of ( 11 | 0 || 1 ) and ( 0 || 0 || 0 ), while ( 2 || 2 || 0 ) 
is just the arithmetical average of (1||1||0) and (0||0||0). This point 
is apparent in Eig. p[ if we notice that in each group the three curves 
intersect at one point. Second, applying the Wigner-Eckart theorem 
for the cubic group, we can write several linear identities consisting 
of the lattice matrix elements. They concern not only the compo¬ 
nents shown in Eig. ^ but also the one that vanish as a — 7 > 0. 

In the upper panel of Eig. ^ most of the components show oscil¬ 
lations and have more than one extremum in this region. Eor exam¬ 
ple, ( 0 || 0 || 0 ) reaches a maximum at a = 1.6 fm and two minima at 
a = 1.2 fm and 2.1 fm, respectively, while the ( 1 || 0 || 1 ) only has an 
apparent minimum at a = 1.8 fm. Eor large lattice spacings, the cal¬ 
culated values usually deviate from the continuum ones by 50% to 
100%. If we want to approximate the continuum limit using a sin¬ 
gle component in the lattice calculations, ( 2 || 11 | 1 ) is the best choice. 
The corresponding curve only deviates slightly from the continuum 
value 12.5 fm^ when a >2.0 fm, which is sufficient for most recent 
lattice simulations of nuclei. 

An anomaly occurs at a = 2.1 fm where the calculated matrix el¬ 
ement ( 0 || 0 || 0 ) becomes negative, whereas the continuum limit is 
a definitely positive number (r^). This discrepancy arises because 
we do not calculate the expectation value of the operator with the 
same wave functions on both sides, as was the case for scalar oper¬ 
ators as shown in Eig. [I|. Eor the I = 2 matrix elements in Eig. 
the reduced matrix element ( 0 || 0 || 0 ) is defined to be proportional 
to the expectation value of the quadrupole operator r^l 20 - On the 
lattice the angular part of the quadrupole operator can not be sepa¬ 
rated completely, thus the insertion of the spherical harmonic T 20 is 
not fully canceled by the C-G coefficients included in Eq. Jl^). The 
remaining lattice artifacts may become as large as the magnitude 
of the observable itself. Therefore, it is dangerous to use randomly 
selected matrix elements calculated at large a as approximations of 
the continuum results, otherwise we may even find vanishing results 
for certain lattice spacings, which is unrealistic. 

In spite of the diversity of the behaviors of the various compo¬ 
nents, we can eliminate the anisotropy of the lattice artifacts by 
calculating the isotropic average according to Eq. In the up¬ 
per panel of Eig. the isotropic average as a function of the lat¬ 
tice spacing is represented by the solid line. In this case, the ex¬ 
pression can not be written as a simple multiplet-weighted average 
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Figure 2. (color online). Upper panel: Mean square radii (r^) for the lowest 
2+ states of the ^Be nucleus as a function of a. (o!||/3 ||y) is an abbreviation 
of the reduced lattice matrix element (2a||r^y2/3 II^T) defined in Eq. (]I^, 
which converges to (r^) as a —> 0. The solid line represents the isotropic 
average (2||r^y2 ||2)o dehned in Eq. Lower panel: Mean value {/') for 
the lowest 2+ states of the *Be nucleus as a function of a. (o!||/3||y) is an 
abbreviation of the reduced lattice matrix element ( 2 a||r^y 2/3 II^T); which 
converges to (r^) as a —0. The solid line represents the isotropic average 

(2|k4y2l|2)o. 


because the C-G coefficients are no longer all the same. Similar 
to the isotropic average shown in the upper panel of Fig. the 
isotropic average in Fig. ^ curve slightly bends downward in the re¬ 
gion 1.5 fm < fl < 2.0 fm. For the lattice spacings considered here, 
the range of values obtained are between 9.8 fm^ and 12.5 fm^, ac¬ 
counting for a 20% relative error with respect to the continuum limit. 
The component that is closest to the isotropic average for all lattice 
spacings is (1||2||1). 

Now let us change the radial factor of the inserted operator and 
keep the angular part the same. In the lower panel of Fig. || 
the symbol (a||j3||7) denotes the lattice reduced matrix element 
{2a\\r%p \\2y) which converges to {r'^) as a —>■ 0. The solid line 
represents the isotropic average (2||r^y2||2)o. Again the inserted 
spherical harmonic ¥ 2 ^ is not fully canceled by the C-G coefficients 
on the lattice. The remaining lattice artifacts shift the various com¬ 
ponents from the continuum limit by different amounts for large 
lattice spacings. Nevertheless, comparing the curves denoted by the 
same symbol in the upper and lower panel of Fig. we find that 
their behavior is qualitatively similar. For example, the (0||0||0) 
curves both show a maximum at a = 1.6 fm and a minimum at 2.1 



Figure 3. (color online) Mean value (r‘^) for the lowest 2+ states of the 
®Be nucleus as a function of a. (o!||/l ||y) is an abbreviation of the reduced 
lattice matrix element (2o!||r^y4^ ||2y) defined in Eq. ([1^, which converges 
to (r'^) as a —> 0. The solid line represents the isotropic average (2||r^y4||2)o 
defined in Eq. (^^. 


fm, while the (2||0||2) curves both show a minimum at 2.0 fm. In 
other words, the magnitude of the lattice artifacts may be different 
if the radial part of the inserted operator is changed, but the pattern 
of oscillations and the sign of the deviations are largely determined 
by the angular momenta of the states and the irreducible tensor op¬ 
erator. 

Next, we turn to the irreducible tensor operators with / = 4. In 
Fig. 1^ the symbol (a||j3||7) denotes the lattice reduce matrix ele¬ 
ment ( 2 a||r^y 4 j 31127 ) which converges to {r'^) as a —>■ 0. The solid 
line represents the isotropic average (2||r"’'y4||2)o. Now the number 
of independent components is larger than that in the case of / = 2 
and the situation is different. For example, the (0||0||0) curve is 
much closer to the isotropic average compared to the corresponding 
curve in the lower panel of Fig. while the (2||0||2) curve has a 
pronounced minimum at a = 2.0 fm. The component closest to the 
continuum limit as well as the isotropic average for all lattice spac¬ 
ings considered here is (11|3||2). There also exist a number of linear 
relations among these components, which is a consequence of the 
remaining unbroken symmetry group. There is also a negative com¬ 
ponent (2||0||2) at a = 2.0 fm. We will not discuss all the details 
because most of them are natural extensions of the I = 2 case. 

Finally, let us examine the postulate that the anisotropy result¬ 
ing from the lattice artifacts is removed in the isotropic average, 
Eq. (^^. In the upper panel of Fig. Q we show the comparison be¬ 
tween the isotropic averages (2||r^yo||2)o and (2||r^y2||2)o- The for¬ 
mer one is calculated with a simple multiplet averaging over the five 
fold branches, while the latter one is obtained by a relatively com¬ 
plicated averaging with the C-G coefficients as the weights. Clearly, 
for all the lattice spacings considered here, the difference between 
the two curves is rather small. Especially, the down bending occurs 
for the same lattice spacing and the magnitudes are also similar. Re¬ 
membering that the C-G coefficients are included explicitly in the 
definition of isotropic average Eq. (^^, we conclude that the effect 
of the angular part of the inserted operators is canceled by the C-G 
coefficients in the nominator. As is discussed in the previous sec¬ 
tion, this is a strong evidence for approximate rotational symmetry 
restoration and that we have factorized the radial and angular parts 
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Figure 4. (Color online) Upper panel. Mean square radii (r^) for the low¬ 
est 2+ states of the ®Be nucleus as a function of a. The black and red lines 
represent the isotropic average (2||r^Ko||2)o and ( 2 ||r^y 2 || 2 )o, respectively. 
Lower panel. Expectation value (r^) for the lowest 2+ states of the ^Be nu¬ 
cleus as a function of a. The black, red and blue lines represent the isotropic 
averages (2||r^yol|2)o, ( 2 ||r^y 2 || 2 )o and ( 2 ||r'^y 4 || 2 )o, respectively. 

of the lattice wave function by means of isotropic averaging. 

In the lower panel of Fig. ^ we show the quantities that con¬ 
verge to (r'^) as a —>• 0, including (2||r'^yol|2)o, ( 2 ||r"^l 2 || 2 )o and 
(2||r^}4||2)o. The three curves coincide with each other even for 
a >2.0 fm, which means that the rotational symmetry is restored to 
a large extent after the isotropic averaging. In particular, the differ¬ 
ence between the / = 2 and I = 4 results is almost zero for all the 
lattice spacings. 

IV. SUMMARY 

In lattice calculations, discretization errors can break rotational 
symmetry. The degeneracy of bound state multiplets with the same 


angular momentum is broken on the lattice. The resulting wave 
functions are classified according to the irreducible representations 
of the cubic group instead of the full SO(3) rotational group. For 
most of the observables represented by the irreducible tensor opera¬ 
tors, the relations between the various components become compli¬ 
cated. In this paper, we used an a cluster model to investigate the 
lattice matrix elements of such operators. We have shown that the 
qualitative behaviors of the various matrix elements versus lattice 
spacing is mainly determined by the angular momenta of the states 
and operators. The matrix elements of different operators with the 
same angular momentum show a similar behavior as functions of 
the lattice spacing. 

We have also defined an isotropic average in which the vari¬ 
ous components of the matrix element are linearly combined. The 
weight for a component is proportional to the corresponding C-G 
coefficients with the same quantum numbers. We have shown that 
such a calculation is equivalent to averaging over all possible lattice 
orientations. In such a calculation, we eliminate to a good approxi¬ 
mation the anisotropy of the lattice artifacts. This point is illustrated 
by numerical calculations for the ®Be nucleus with the a cluster 
model. 

We have calculated the isotropic averages as functions of the lat¬ 
tice spacing for irreducible tensor operators with different angular 
momenta and same radial factors and find good agreement with con¬ 
tinuum limit values. Comparing the results calculated with different 
lattice spacings, we found that there is a down bending of (r^) and 
(r^) in the region 1.7 fm < a < 2.0 fm. For a < 1.7 fm, the lattice 
artifacts are very small. 

Although all the conclusions in this paper are obtained with a 
simple a cluster model, the results can be applied immediately to 
ab initio lattice simulations as well. For example, when calculating 
the transition amplitude between low-energy excited states and the 
ground state of a nucleus, we can improve the results by calculat¬ 
ing the isotropic average to eliminate rotational symmetry breaking 
effects on the lattice. 
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